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Abstract 

We report the first experimental observation of extreme muhistabihty in a controhed laboratory 
investigation. Extreme multistability arises when infinitely many attractors coexist for the same 
set of system parameters. The behavior was predicted earlier on theoretical grounds, supported by 
numerical studies of models of two coupled identical or nearly identical systems. We construct and 
couple two analog circuits based on the Lorenz system and demonstrate the occurrence of extreme 
multistability through the controlled switching of the coupled system to different attractor states 
purely through a change in initial conditions for a fixed set of system parameters. Numerical 
studies of the coupled Lorenz circuits are in agreement with our experimental findings. 
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Multistability is a common occurrence in many nonlinear dynamical systems, correspond- 
ing to the coexistence of more than one stable attractor for the same set of system parameters 
[1]. A large number of theoretical and experimental studies have explored this phenomenon 
in a variety of physical [2-5], chemical [6, 7] and biological [8, 9] systems. A curious and 
novel manifestation of this phenomenon arises when a system can have an infinite number 
of coexisting attractors, where each attractor is associated with a particular set of initial 
conditions [10-12]. This extreme multistability has been found in coupled identical or nearly 
identical chaotic systems, such as the Lorenz system [13] and the three-variable autocata- 
lator model [14]. The coupled subsystems exhibit generalized synchronization and, in the 
case of an experimental system with unavoidable variability in conditions, a sensitivity to 
initial conditions with a corresponding 'uncertainty' in the destination dynamics. Extreme 
multistability was first demonstrated in a system of two coupled identical Lorenz oscillators 
by Sun et al. [10] and has more recently been studied in the three-variable autocatalator 
model by Ngonghala et al. [12]. 

In both studies, a special coupling was applied between two three-variable chaotic 
systems to form six-variable coupled systems. Numerical simulations of the coupled systems 
were carried out for a fixed set of system parameters and only the initial conditions were 
changed. The systems were found to evolve to different attractor states (fixed points, limit 
cycles, chaotic states) purely through changes in initial conditions, typically with a change 
in the initial condition of just one of the state variables. This behavior was explained 
theoretically [10, 12] by analyzing the dynamics of a reduced set of equations governing the 
differences ("errors") of the corresponding state variables of the three- variable subsystems. 
The reduced system yields a constant c, a conserved quantity that attains a particular value 
following the decay of transients. Different initial conditions gave rise to new dynamical 
behavior and a corresponding different value of c. The existence of this conserved quantity 
leads to a 'slicing' or foliation of the state space into manifolds corresponding to the 
different values of c. The conserved quantity can be used to reduce the six-variable coupled 
system to a three-variable system, in which c serves as a bifurcation parameter. Although 
the full coupled system has a conserved quantity, the asymptotic states to which the system 
evolves are true attractors, since there are an infinite number of initial conditions that yield 
any particular value of the conserved quantity c. The set comprising these initial conditions 
constitutes the basin of attraction for a specific attractor, and trajectories emanating 
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from such initial conditions will converge to the attractor. Extreme multistability might 
have important consequences in the reproducibility of certain experimental systems. For 
example, some chemical reactions, such as the chlorite-thiosulfate reaction [15] and the 
chlorite-iodide reaction [16] consistently exhibit irreproducibility: despite great care to 
ensure reproducibility, these reactions show a random long-term behavior for the same 
set of experimental conditions. The cause of the irreproducibility is not known; however, 
extreme multistability offers a possible mechanism for the behavior. Model calculations by 
Wang et al. [11] demonstrated a type of extreme multistability arising from the specific 
features of a chemical reaction. However, to the best of our knowledge, there has yet been 
no direct experimental verification of this new type of dynamical behavior in a controlled 
laboratory investigation. In this Letter, we report experimental observations of a coupled 
electronic circuit system that displays extreme multistability. 



Our experiments are carried out on an analog circuit system based on the original model 
investigated by Sun et al. [10], consisting of a set of coupled Lorenz equations, namely. 



Xi = aiYr - X2) (la) 

Yi = rXi -Yi- XiZi (lb) 

Zi =XiYi- bZi (Ic) 

X2 = a(r2 - X2) (Id) 

Y2 = rXi -Y2- (le) 

Z2 =XiY2- bZ2 (If) 

where a, r and b are constants and (Xi, Yi, Zi) and {X2, Y2, Z2) are the state variables of 
the two identical subsystems. For the circuit implementation of the above system, we have 
scaled the state variables as 

_ -^1,2 _ ^1,2 Z12 

^1,2 — /— - ; yi,2 — r— \ Zi2 



/3r ' V3r ' 3r 
in order to restrict the output signal voltage range to within ±10 volts and to avoid saturation 
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FIG. 1: Circuit diagram of two coupled Lorenz oscillator systems 
of the circuit. This leads to the following system: 

(2a) 
(2b) 

(2c) 
(2d) 
(2e) 

(2f) 

The detailed circuit diagram of the two coupled systems is shown in Fig. (1), in which we 
have adopted the basic circuit design for a single Lorenz oscillator recently suggested by N. 
J. Corron [17]. A regulated power supply of ±15^ energizes the circuit, and the system 
parameters r, a, and b for the individual Lorenz circuits are controlled with three circuit 
resistors, namely, 

lOkn lOkQ loom lOOkQ , lOOkQ lOOkQ 

r = = ; cr = = ; o = = 

3i?3 3i?i3 i?7 i?i7 R2 R12 

As a benchmark exercise, each Lorenz oscillator was separately tested by varying the system 
parameters to obtain its various attractor states, ranging from periodic states to chaotic 
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FIG. 2: Timer relay circuit for setting initial conditions of some the state variables to zero. 

dynamics, and the behavior was then compared to numerical simulations of the Lorenz 
equations. Major care was taken to ensure that the two oscillator systems were as nearly 
'identical' as possible within practical limits. This entailed careful weaning of all the com- 
ponent elements (resistors, capacitors) to match their values as closely as possible and the 
removal of any intrinsic drifts or biases within the operational amplifiers and multipliers 
[18]. The system parameters for both of the oscillators were then fixed at cr = 10, r = 400, 
6 = 8/3, and they were coupled to each other in the manner shown in Fig. (1). 

To change the initial conditions of the dynamics of the circuit, we have employed a 
strategy of imposing short voltage pulses on selective nodes of the operational amplifiers as 
well as shorting relevant capacitors of the circuit initially to set the values of some of the 
state variables to zero. To implement this combined strategy in a controlled manner we have 
developed and attached an additional system of relay circuits to the coupled Lorenz circuits 
(indicated by the boxed hatched portions in Fig. (1)). The details of the relay circuit are 
shown in Fig. (2). When this circuit is energized, the timer portion of the circuit produces 
a high output for 1 second and the relay is turned on through the relay driver circuit. Due 
to this, the y2 and ?/2m signals get shorted in the coupled circuit and hence y2 is set to 
initially. Similarly zi and Z2 can also be set to 0. After 1 second the output of the timer 
circuit falls to a low value and the relay gets switched off. Due to this, the y2 and |/2m signal 
paths open up and the circuit can run with the applied initial conditions. Note that we can 
also allow y2 (and or yi) to have a finite voltage initially by not shorting it. In this case 
the initial value is given by the output of the multiplier and op-amp and can be determined 
exactly from the circuit parameters. The initial conditions of xi and X2 can be changed 
through application of short voltage pulses using the VI and VA voltage sources as shown 
in the circuit diagram. 
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FIG. 3: The right panels show snapshots of experimental attractor states of the coupled electronic 
circuit implementation of the model system, Eqs. (2). The oscilloscope images show Z2 vs xi phase 
space plots for different initial values. The left panels show corresponding plots from numerical 
simulations of Eqs. (2) with similar initial conditions. 

Using the above strategy we have run the coupled circuit for a number of initial conditions 
without changing the circuit parameters. Some typical results in the form of oscilloscope 
images of phase plots of Z2 vs Xi are shown in the right-hand column of Fig. (3) indicating two 
different periodic attractor states and a chaotic state. The uppermost image corresponds 
to initial conditions of (xi = +2V,yi = Zi = X2 = y2 = Z2 = 0), the middle image to 
(xi = —2V, yi = zi = X2 = y2 = Z2 = 0) and the lower image has (xi = — 2V, X2 = —IV, yi = 
z\ = Z2 = 0,y2 = 14.2V). The plots in the left-hand column show corresponding numerical 
simulation results of the model Eqs. (2) using the same initial conditions. In Fig. (4), we 
show the oscilloscope wave forms of the variables Xi, X2 and the time variation of their 
difference (xi — X2). In all three cases, the differences attain constant values of different 
magnitudes corresponding to the three distinct attractor states, which is a signature of 
extreme mult ist ability. The difference (xi — X2) represents the conserved quantity c in the 
system. 

The coupled Lorenz circuit studied here, as well as other coupled systems that exhibit 
extreme multistability, exhibit properties of both dissipative and conservative dynamics. 
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FIG. 4: Experimental waveforms of xi, X2 and X1 — X2, where the left and central panels correspond 
to the periodic attractor states and the right panel to the chaotic attractor state as shown in Fig. (3). 

The circuit is characterized by infinitely many attractors, each associated with a particular 
value of the conserved quantity c = xi — X2, where the basin of attraction is made up of 
all sets of initial conditions that evolve to asymptotically yield the particular value of c 
associated with the attractor. Similarly, the attractor is stable to perturbations in which 
the asymptotic value of the conserved quantity is maintained. As pointed out in Ref. [12], 
because the manifold of initial conditions associated with any attractor is arbitrarily close to 
manifolds of initial conditions associated with other attractors, all attractors in our system 
are weak attractors in the Milnor sense [19]. 

The conserved quantity gives rise to a direction of neutral stability for the stationary 
states as well as the orbits (in addition to that associated with the direction along the or- 
bit). Hence, if the system is in a particular periodic orbit, say period-2, a perturbation that 
does not satisfy the condition of the conserved quantity will give rise to the evolution of the 
system to a new attractor. The new attractor may differ only quantitatively; for example, a 
small perturbation might shift the period-2 dynamics to a new period-2 dynamics that differs 
in amplitude. However, larger perturbations give rise to the evolution of the system to quali- 
tatively new attractors, such as a period-4 or period-8 attractor. This characterization of the 
effects of perturbations also applies to the effects of different initial conditions. Because the 
coupled system exhibits period-doubling chaos, perturbations or different initial conditions 
permit the sampling of any of an infinite number of qualitatively different attractors. Even 
if the coupled system did not display chaotic dynamics, the same mechanism would give rise 
to an infinite number of quantitatively different attractors. As yet, however, we have been 
unable to find the phenomenon of extreme multistability in non-chaotic coupled systems. 

We have described the first experimental demonstration of extreme multistability. Our 
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experiments have been based on an electronic circuit that models two coupled Lorenz at- 
tractors. A restrictive feature of this new type of dynamical system is the requirement that 
the coupled subsystems be identical or nearly identical [12]. It is this feature that makes 
electronic circuit models especially attractive, since, as we have demonstrated, two chaotic 
circuits can be sufficiently matched to give rise to extreme multistability when appropriately 
coupled. It is likely that extreme multistability will be a rarity in most physical, chemical and 
biological systems; however, the combined conservative and dissipative features give rise to 
dynamics that might find technological uses, such as the ability to easily select qualitatively 
different dynamical states from an infinite number of possibilities. 
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